Beyond the Shell Model: The Canonical Nuclear Many-Body Problem as an Effective 
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We describe a strategy for attacking the canonical nuclear structure problem — bound-state prop- 
erties of a system of point nucleons interacting via a two-body potential — which involves an expan- 
sion in the number of particles scattering at high momenta, but is otherwise exact. The required 
self-consistent solutions of the Bloch-Horowitz equation for effective interactions and operators are 
obtained by an efficient Green's function method based on the Lanczos algorithm. We carry out 
this program for the simplest nuclei, d and 3 He, to contrast a rigorous effective theory with the shell 
model, thereby illustrating several of the uncontrolled approximations in the latter. 



In this letter we argue that it may be possible to move 
beyond the nuclear shell model (SM) to a more rigor- 
ous treatment of the canonical nuclear structure prob- 
lem of A nonrelativistic point nucleons interacting via a 
two-body potential. Our optimism is inspired by several 
recent developments. One is the success of the Argonne 
group's efforts JlJ] to predict the properties of light nuclei 
in effectively exact Green's function Monte Carlo calcu- 
lations, using an NN potential carefully fit to scattering 
data and augmented by weaker three-body forces. This 
suggests that SM failures have their origin in an incom- 
plete treatment of the many-body physics, rather than 
in the starting Hamiltonian. A second is the success of 
effective field theory (EFT) treatments of the two- and 
three-body problems. This work not only provides some 
insight into why such a starting Hamiltonian is reason- 
able, but has made the community more aware of the un- 
controlled approximations implicit in the SM and other 
approaches. 

The SM's strength is its explicit representation of ~ 
60% of the wave function that resides at long wave- 
lengths: the A-body correlations important to collec- 
tive modes are addressed by direct diagonalization. A 
third development is the remarkable recent advances in 
such SM technology, including Lanczos-based methods 
M , treatments of light nuclei involving many shells @, 
and Monte Carlo sampling algorithms The dimen- 

sions of tractable SM spaces have risen by several orders 
of magnitude in the past few years. 

Such diagonalizations in a long-wavelength "included 
space" could be an important piece of a rigorous effective 
theory (ET) of nuclear structure in which the Hamilto- 
nian operating in an infinite Hilbert space 
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is replaced by an H & '' operating in a finite SM space. 
The effects of high momentum components appear as ef- 
fective contributions to the Hamiltonian and operators. 
The hope in nuclear physics, inspired by Brueckner's 



treatment of nuclear matter, is that the "excluded space" 
integration might be carried out as a rapidly converging 
series in the number of nucleons scattering at one time in 
high momentum states. In this way the effective theory 
might prove far more tractable than the original A-body 
problem in a infinite Hilbert space. 

Despite some text book motivations, the SM is a model 
rather than an ET: 

1) The functional form of the SM effective interaction, 
(a[3\H e ff\j5), is correct only in lowest order and only 
if the calculation is restricted to a single shell In a 
faithful ET three-, four-, and higher-body operators are 
successively added, and the matrix elements generally 
carry, in addition to single-particle quantum numbers, 
an index specifying the number of quanta carried by the 
remaining spectator nucleons. 

2) Typically H e ^ lacks the symmetries of the original 
bare H, e.g., translational invariance and hermiticity 
(though the latter is often enforced by hand). 

3) SM wave functions are orthogonal and normed to 
unity. In ET the effective wave functions are naturally 
defined as the restrictions of the true wave functions 

to the model space, so that the norms are less than unity 
and orthogonality is lost. 

4) Shell model interactions frequently depend on ficti- 
tious parameters such as "starting energies," introduced 
to adjust the energy denominator in the two-body G- 
matrix or to account for intermediate-state average ener- 
gies when the two-body G-matrix is iterated to produce 
some approximation to a higher order H e ''. 

5) Perhaps most serious, the important issue of effective 
operators is almost never addressed in a meaningful way. 
In most cases SM H e ^s are derived phenomenologically, 
so that there is no diagrammatic basis for generating the 
effective operator. Thus empirical operator renormaliza- 
tions must also be introduced, obscuring the underlying 
physics and undercutting the SM as a predictive tool. 

In this letter we describe the first steps in an effort to 
assess the feasibility of an exact ET "SM"-like theory. 
The approach is sketched in Fig. 1. The Hilbert space 
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FIG. 1. Cluster expansion of the effective interaction. 

is divided into a long-wavelength "SM" space, defined by 
some energy scale Asm, and a high- momentum space. 
One can truncate the latter at some scale A^ ~ 3 GeV, 
characteristic of the cores of realistic potentials, as above 
this energy, excitations make a negligible contribution. 
All correlations within the "SM" space arc included, but 
the high-momentum correlations in the excluded space 
are limited to n-body, where n is the cluster size. Thus 
Figs, lb and lc give the lowest and next-to-lowest ap- 
proximations to H e ff. Note that one can view Fig. lc 
as containing Fig. lb, a density-dependent two-body cor- 
rection to Fig. lb, and true 3-body terms. The pattern 
continues as n is increased: true n-body terms are in- 
troduced and all lower-order results are corrected to one 
higher order in the density p. 

If the "SM" space is defined as all harmonic oscillator 
Slater determinants with E < As m , H e ->* becomes trans- 
lationally invariant and the ladder sums can be carried 
out in relative coordinates, a considerable simplification. 
The projection operator onto the high-momentum space 
Q thus depends on Asm and the oscillator parameter 6, 
where the latter can be chosen to optimize the conver- 
gence in Aoo 0. 

The resulting Bloch-Horowitz (BH) equation [pj is then 



H eff = H + H- 
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H eit \* S m) = E\* S m) |*sm) = (1 - Q)|* ) (2) 

where |\&) is the exact wave function and H\^>) = E\^f). 
These equations must be solved self-consistently because 



H e ff depends on the unknown eigenvalue E. The har- 
monic oscillator appears only implicitly through Q in 
distinguishing the long-wavelength "SM" space from the 
remainder of the Hilbert space. 

There is an extensive literature on this and similar 
equations ||Jl(|. Frequently H is divided into an unper- 
turbed Hq and an perturbation H — Hq, but well-known 
pathologies due to intruder states can affect the result- 
ing perturbation expansion |ll]] . The approach explored 
here is nonperturbative and uses the Lanczos algorithm 
to sum the n-body ladders. 

The Lanczos algorithm recursively maps a hermitian 
operator H of dimension N into tridiagonal form 

H\vi) =ai|«i) + /3ih) 

H\v 2 ) = /8i|«i)+aa|«2)+/%K) 

H\v 3 )= f3 2 \v 2 ) + a 3 \v 3 ) + /3 3 \v 4 } etc. (3) 

If this process is truncated after n <C N steps, the re- 
sulting matrix contains the information needed to recon- 
struct the exact 2n-l lowest moments of \v±) over the 
eigenspectrum. One of the applications of this algorithm 
is in constructing fully interacting Green's functions |l2j ] 
function of E, 
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where the gi{E) are continued fractions that depend on 
cti,(3i and where E appears only as a parameter. 

Thus a simple procedure can be followed to solve the 
BH equation: 

• For each relative-coordinate vector in the SM space 
1 7), form the excluded-space vector \vi) = QH\j) and 
the corresponding Lanczos matrix for the operator QH . 
Retaining the resulting coefficients ccj,/3j for later use, 
construct the Green's function for some initial guess 
for E and then the dot product with (j'\H to find 
( 7 '|iF//(£)|7>. 

• Perform the "SM" calculation to find the desired eigen- 
value E' which, in general, will be different from the guess 
E. Using the stored OLufii, recalculate the Green's func- 
tion for E' and H eif (E') then redo the "SM" calcula- 
tion. The process is repeated until the energy is fully 
converged. 

• Then proceed to the next desired bound state and re- 
peat the process. Note that it is not necessary to re- 
peat the H e ff calculation. The eigenvalue taken from 
the "SM" calculation is, of course, that of the next de- 
sired state, yielding a distinct H e *f for each eigenvalue. 

The attractiveness of this approach is that the effec- 
tive interactions part of the procedure, which is relatively 
time consuming as it requires one to perform a large-basis 
Lanczos calculation for each relative-coordinate starting 
vector in the "SM" space, is performed only once. While 
the SM calculation must be repeated in the iterations, the 
convergence is rapid (6-8 cycles, typically). As modern 
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workstations can manage calculations of dimension 10 6 
in about 30 minutes, the overall procedure is certainly 
practical. 

The technical aspects of this approach are described 
elsewhere [7|JT^|. Here we focus on the results for the 
simplest nuclei, d and 3 H, carrying the above process to 
completion (two- and three-body ladders, respectively). 

The binding energies and operator matrix elements for 
simple systems like 3 He can, of course, be calculated ex- 
actly by other methods. The point of our work is not 
to offer an alternative to these techniques for these nu- 
clei, but rather to illustrate the conceptual differences 
between a faithful ET and the SM. We performed d and 
3 He ET calculations for a series of "SM" spaces (2, 4, 
6, and 8 7kS), in each case using the Lanczos Green's 
function algorithm to evaluate the two- and three-body 
ladders (100 Lanczos iterations are more than sufficient) 
and iterating the "SM" calculation until the results are 
fully converged. The deuteron calculation is rather triv- 
ial; for Aqo ~ 50 the 3 He calculation involves a dense 
matrix of ~ 10 4 , still quite modest by current SM stan- 
dards. (The Hamiltonian matrix is dense because relative 
Jacobi coordinates are used, rather than the m-scheme, 
together with standard Talmi-Brody-Moshinsky methods 
|^3|,0.) A first test of this procedure is the stability of 
the resulting energy eigenvalue: the results for the four 
chosen "SM" spaces agreed to four places, -2.224 MeV 
(using \[2b =1.6/ and Aoo = 140). The exact result is 
-2.2246 MeV. The 3 He results were similarly very stable. 

A more interesting test is the evolution of the wave 
function as the "SM" space is enlarged. Table | gives 
results for 3 He. (The procedure for calculating the wave 
function normalization is discussed below.) Unlike typ- 
ical SM calculations, the amplitudes agree over over- 
lapping pieces of the Hilbert spaces. As one proceeds 
through 2fiu!, 4huj, 6%w, ... calculations, the ET wave 
function evolves only by adding new components in the 
expanded space. Consequently, as Table I shows, the 
wave function norm grows. 

This evolution will not arise in the standard SM be- 
cause the wave function normalization is set to unity re- 
gardless of the model space. It will also not arise for a 
second reason, illustrated in Table |l[ The three-body 
3 He matrix elements of H e * ' are crucially dependent on 
the model space: a typical matrix {a\H e *' changes 
very rapidly under modest expansions of the model space, 
e.g., from 2Tiw to ATiuj. Yet it is common practice in the 
SM to expand calculations by simply adding to an ex- 
isting SM Hamiltonian new interactions that will mix 
in additional shells. We suspect the behavior found for 
3 He is generic in ET calculations: it arises because a 
substantial fraction of the wave function lies near but 
outside the model space (e.g., see Table ||). An expan- 
sion of the model space changes the energy denominators 
for coupling to some of these configurations, and moves 
other nearby configurations from the excluded space to 



TABLE I. ET results for the 3 He ground state wave func- 
tion calculated with the Argonne v 18 potential. Selected ba- 
sis states are designated somewhat schematically as | N,a), 
where N is the total number of oscillator quanta and a is an 
index representing all other quantum numbers. 
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amplitude 


Qhuo 






Qfno 


8 huj 


exact 


(31.1%) 


(57.4%) 


(70.0%) 


(79.8%) 


(85.5%) 


(100%) 


0,1) 


0.55791 


0.55791 


0.55791 


0.55795 


0.55791 


0.55793 


|2,1) 


0.00000 


0.04631 


0.04613 


0.04618 


0.04622 


0.04631 


1 2,2) 


0.00000 


-0.48255 


-0.48237 


-0.48243 


-0.48243 


-0.48257 


1 2,3) 


0.00000 


0.00729 


0.00731 


0.00730 


0.00729 


0.00729 


I 4,1) 


0.00000 


0.00000 


-0.02040 


-0.02042 


-0.02043 


-0.02047 


1 4, 2) 


0.00000 


0.00000 


0.11267 


0.11274 


0.11275 


0.11289 


1 4, 3) 


0.00000 


0.00000 


-0.04191 


-0.04199 


-0.04208 


-0.04228 



TABLE II. Selected BH 3-body effective interaction matrix 
elements for 3 He, in MeV, illustrating the strong dependence 
on the "SM" space. 





2hcu 




6hui 


8 huj 


(0, 1 | H"~ n | 2, 1) 


-4.874 


-3.165 


-0.449 


1.279 


(0, 1 | H eti | 2, 5) 


-0.897 


-1.590 


-1.893 


-2.208 


(2,1 H ets | 2,2) 


6.548 


-2.534 


-4.144 


-5.060 



the model space. Naively, relative changes in effective 
interaction matrix elements of unity are expected. 

Now we turn to the question of operators. The stan- 
dard procedure in the SM is to calculate nuclear form fac- 
tors with bare operators, or perhaps with bare operators 
renormalized according to effective charges determined 
phenomenologically at q 2 = 0, using SM wave functions 
normed to 1. As we now have a series of exact effective 
interactions corresponding to different model spaces, we 
can test the validity of this approach. The results for 
the elastic magnetic form factors are given in Fig. 2. 
Even though each "SM" H e ^ is, in a sense, perfect, the 
results for bare operators are widely divergent at even 
modest momentum transfers of ~ 2/f. This is not sur- 
prising: if an operator transfers a momentum q > 2kp 
to the nucleus, where Uf is the Fermi momentum, the 
resulting amplitude should reside primarily outside the 
"SM" space, where it contributes only through effective 
pieces of the operator. 

Clearly the effective interaction and effective operator 
have to be treated consistently and on the same footing. 
The bare operator O must be replaced by 



O eff = (1 + HQ- 



1 



;)0(1 



1 



;QH) (5) 



Ef - HQ v Ei-QH 
and must be evaluated between "SM" wave functions 
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FIG. 2. The magnetic elastic form factors for the deuteron 
(top) and 3 He (bottom) calculated with the exact H e ^ , SM 
wave functions normalized to unity, and a bare operator are 
compared to the exact result (solid line). When effective op- 
erators and the proper wave function normalizations are used, 
all results become identical to the solid line. 



normed according to 
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(and similarly for |\E i y SM )). These expressions can be 
evaluated with the Lanczos Green's function methods de- 
scribed earlier. When this is done, all of the effective 
calculations, regardless of the choice of the model space, 
yield the same result, given by the solid lines in Fig. 2. 

It is likely that many persistent problems in nuclear 
physics — such as the renormalization of in P decay 
— are due to naive treatments of operators. It should be 
apparent from the above example that no amount of work 
on H e M will help with this problem. What is necessary 
is a diagrammatic basis for generating H e ff that can be 
applied in exactly the same way to evolving O e ^ . From 
this perspective, phenomenological derivations of H e *' 
by fitting binding energies and other static properties of 
nuclei are not terribly helpful, unless one intends to si- 
multaneously find phenomenological renormalizations for 
each desired operator in each q 2 range of interest. 

There are several remarks we would like to make in 
conclusion. First, the demonstration that one can ef- 
ficiently solve the three-body problem as an exact ET 



in a SM-like model space is already quite significant: 
this means it is relatively straightforward to execute a 
faithful BH treatment of heavier nuclei through order 
p in both the effective interaction and effective opera- 
tors. The numerical effort is comparable. Second, our 
results indicate the many of the approximations in the 
SM are uncontrolled, so that a more faithful ET is badly 
needed. Third, further progress can be made numeri- 
cally: we expect to gain a factor of 50 when the large- 
memory capabilities of our workstation are properly ex- 
ploited. Finally, there is the strong possibility of theoret- 
ical improvements. For example, the numerical conver- 
gence of our binding energies for large Aqo is exponential, 
exp(— aA^,). We argue in J?j that the origin of this scal- 
ing is the contraction of ladders at very high excitation 
energies to local operators, which can be summed ana- 
lytically in a local density expansion. If this conjecture 
proves true, we can bring Aoo down to a much lower scale, 
with no loss of precision. Thus much remains to be done. 
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